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Abstract. - We examine the practical feasibility of the experimental realization of the so-called 
entangled Bose-Einstein condensation (BEC), occurring in an entangled state of two atoms of 
different species. We demonstrate that if the energy gap remains vanishing, the entangled BEC 
persists as the ground state of the concerned model in a wide parameter regime. We establish the 
experimental accessibility of the isotropic point of the effective parameters, in which the entangled 
BEC is the exact ground state, as well as the consistency with the generalized Gross-Pitaevskii 
equations. The transition temperature is estimated. Possible experimental implementations are 
discussed in detail. 
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Recently, a novel type of Bose-Einstein condensation 
(BEC), characterized by an entangled order parameter, 
herein simply called entangled BEC, was found to be the 
ground state of a mixture of two species of atoms with 
spins [1]. The entangled BEC, which occurs in an in- 
terspecies two-particle entangled state, is interesting in 
several aspects. First, it is not a mean field state of sin- 
gle atoms, unlike the usual multicomponent BEC, whose 
ground state is simply a direct product of wave func- 
. tions of all the atoms, with each component described 
■ by a same wave function. Second, the interspecies two- 
particle entangled state bears some similarity to Cooper 
' pairing of fermions, but it is not a weakly bound state. In 
other words, the entangled order parameter is nonlocal. 
Third, this entangled condensate may serve as a source of 
Einstein-Podolsky- Rosen pairs, the key components in a 
lot of quantum information processes. The nonlocal na- 
ture of the order parameter makes it convenient to be used 
for this purpose, if experimentally realized. Fourth, as 
elaborated below, the entangled BEC is a new kind of 
fragmented condensation, i.e., with macroscopic occupa- 
tion of more than one single-particle state [2], which has 
attracted a lot of attention in recent years [3,4]. 

The entangled BEC was previously established as the 



('''Email: yushi@fudan.cdu.cn 



ground state of the concerned model, when the effective 
parameters are so chosen that the total pseudospin of the 
system is conserved, in consistency with the common wis- 
dom that fragmentation is due to spin symmetry [4]. An 
important question is how this ground state can survive in 
a wide parameter regime, where the Hamiltonian breaks or 
lowers the symmetry. This is crucial for experimental re- 
alization. There have been some numerical investigations 
on this issue demonstrating the persistence of interspecies 
entanglement [1]. But there has been a lack of analytical 
discussions. Moreover, it has not been clear enough how 
to actually implement the entangled BEC in experiments. 

In this Letter, we investigate several related issues 
around the experimental realization of the entangled BEC. 
First we consider the concerned model with a deviation 
from the symmetric parameter point, where the entangled 
BEC is the exact ground state. Treating the deviation as a 
symmetry breaking perturbation, we find that for a finite- 
volume gas, if the energy gap tends to vanish, then the 
ground state approaches the unperturbed entangled BEC. 
This establishes a significant parameter regime for the en- 
tangled BEC. Afterwards, we examine the consistency be- 
tween the isotropic point of the effective parameters and 
the generalized Gross-Pitaevskii equations of the orbital 
wave functions, which are in turn related to those effec- 
tive parameters. After a brief estimation of the transition 
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temperature, we discuss in detail how to experimentally 
implement the entangled BEC. In view of recent related 
experimental progress, we suggest that it is feasible to ex- 
perimentally realize the entangled BEC in an optical trap. 

The model. - Consider a dilute gas of two 
species (or isotopes) of bosonic atoms in a trap [1]. 
Each atom possesses an internal degree of free- 
dom represented as a pseudospin with z-component 
basis states t and J,. Then the many-body Hamil- 
tonian is H = J2i=a,bl^a=tl I ^^^'f't'^i'^'^i'^ + 

J <Prtpl^^^l^^U^fXT3T4i>b<Tsi>aai, where the indepen- 
dent variable r in the field operator ipicr, the single 

(ii) 

particle Hamiltonian hia and the interaction [/cti<t2CT3<t4 
are all omitted for brevity, hi^ = —?^V'i/2mi+Uia, where 
Uia is the trapping potential. The s-wave effective inter- 
action Uil^^^aaaiin - r'j) = {2lTh^i'^^ll^asai/ IJtij)5{'^i " 

= a,b), where = mimj/{mi +mj) is the effective 
mass, S,a^}r.20-3iJ4 is scattering length for an allowed channel 
in which the initial pseudospins of i and j are (74 and (73, 
respectively, while their final pseudospins are ai and <J2, 
respectively. As usual, in ignorance of the depletion, the 
orbital degree of freedom of each atom is constrained in 
the manifold of the single-particle orbital ground states 
Hence ipaa = aa4>aa{va) and ^Abo- = bcr4>ba{n), 
where and are annihilation operators. Under the 
conservation of z component of the total pseudospin in 
each scattering, the many-body Hamiltonian is simplified 
to 

+ y,K^aa'NaaNba' + ^{a\a^h\b^ + a\a^b\b^), (1) 

where Ni^r is the number of atoms of species i with pseu- 
dospin a, Ni = Ni-f+Nii is conserved, fia = ei^-K^^a /2, 



-I' 



/2mi + Uu 



(2) 



is the single particle energy of an atom of species i and 
pseudospin a. Note that there is no a priori requirement 



on the relation between 



and 



K, and K^^) are ef- 



fective parameters proportional to the corresponding scat- 
tering lengths, and are defined in the following way. First, 
for the scattering in which an i-atom fiips from ua to a\ 
while a j-atom fiips from 0-3 to (T2, define 



2nh\i% 



d^r. (3) 



Then in accordance with the convention of Leggett [5], 



we use the shorthands K, 



(ii) 



J^a-cra-a-j ancl 



' ' i\ = 2K^^^g for a ^ a for intraspccies scatter- 
ing, while K^^} = k'^^^}^,^ for interspecies scattering. 



For interspecies pseudospin exchange scattering, we de- 



note = 2K 



'TiTi ~ ^"^iUr '^^^ -^ast term of the Hamil- 
tonian (1) represents pseudospin-exchange scattering be- 
tween two atoms of different species. Without loss of gen- 
erality, suppose Na > Nb- 

Ground State cind fragmentation. - Using the total 
pseudospins of the two species, 

Sa = ^ ] (4r^(T(7' O'er' J S(, = ^ ] ^gScrcr' ^tr' , 



and subtracting a constant, the Hamiltonian can be 
rewritten as 

n = no + ni, (4) 

with 

Ho = KeSa ■ Sb, 
= {Jz — Ke)SazSbz + BaSaz + BbSbz + CaS^z + ^bSlz, 

where 



7 _ I<'(af') 
Uz — 



T^{ab) _ ry-iab) _ T^{ab) 

'^ii -^n, -^iT , ' 



Ba = fa] - fa[ + 



■II 



Bb 



AT-/H + f(<'-^U 



T^(ab) _ T^{ab) _ j^(ab)^ 

-^n -^iT -^ii 

(bb)^ , 



2 (-^TT 



Cb = UKf}'^+K 



T^(ab) _ T^(ab) 

(aa) 

Ti, 



T^(ab)^ 



+ Kf^"-^ - K 

bb) _ r.(bb) _ j^ibbh 

2V-"-TT "^-""U Ti iT ''■ 



bi) 



K. 



(bb) 



K 



{bb). 



Ho is symmetric under SU{2S+ 1), where S is the total 
spin quantum number. The eigenstates of Ho arc \S, 5'^)'s, 
where S = Sa - Sb, - ■ ■ ,Sa + Sb, -5, • • • , 5 is the 

^-component of the total pseudospin, Sa = Na/2, Sb = 
Nb/2. Obviously, except the ground state in the case of 
Na = Nb, the eigenstates are degenerate with respect to 

Sz. 

The ground state of Ho, in the sector of Sz, is 

\GsJ=A{a\r{alr^{a\b\-a\b\f^\0), (5) 

for which S = Sa — Sb = {Na — Nb) /2, A is the normaliza- 
tion constant, = Na/2-Nb/2 + Sz, = Na/2-Nb/2- 
Sz, Sz = -{Na-Nb)/2,-{Na-Nb)/2 + l, • • • , {Na-Nb)/2 
is a conserved quantity. 

The ground state \Gs^) is a special kind of fragmented 
condensate. Because a and b are two different species of 
atoms, one should define a one-particle reduced density 
matrix for each species respectively. For species i, it is 



p,(r,r')^(G|it(r>,(r)|G) 



where 



ipaiv) = a.^4>a\(j)\ T) + ai0ai(r)| i), 
V-feW = ^T<^feT(r)l ^)+bl<Pbl{v)\ i). 
It is evaluated that 
p„(r,r') = {Nal2 + Sz)ra]{v')<i>a]{v) 

+ {Nal2-SMl^{v')<j,al{v), 

Pb{v,v') = Nbl2ct>l^(v')4>b]{v) + Nb/24>l^{v')ct>bi{v). 
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Therefore, all atoms of species a form a fragmented 
condensate while all atoms of species h form another 
one. There is one-particle off-diagonal long-range order 
(ODLRO) within each species. 

Note that in our definition of the field operators and 
reduced density matrices, we have included summation 
over the two spin states, i.e. we have used a spinor field 
operator, in considering the coherence between the two 
spin states. The reduced density matrix of each species 
characterize the position order for all the atoms of this 
species. 

The interspecies entanglement leads to nonvanishing 
two-particle reduced density matrix of two particles of dif- 
ferent species, 



i{v',)M^a)Mn)\Gs.), 



which has more than four eigenvalues of order iV^. Hence 
there is ODLRO in the interspecies two-particle density 
matrix. Moreover, there is phase coherence in the inter- 
species pairs. 

The state (5) is the so-called BEG with an entangled 

order parameter, simply called entangled BEG. Note that 
there is ODLRO in the single particle density matrix of 
each species of atoms. 

Persistence of entangled BEC. Away from the 
symmetric point in the parameter space, the symmetry 
breaking perturbation Hi becomes nonvanishing. H still 
conserves 5^, though S is no longer conserved. Let us 
focus on the sector of Sz = 0, in which any eigenstate of 
H can be expanded as 



|*„)=13^n(^)|5,0). 



(6) 



where the summation is over S, running from Smin = 

\Sa — Sb\ to Sa + Sb- 

H|*„) = leads to where N = {Na + Nb)/2, 

E^^^S) = KeS(S+ 1). In deriving this equation, we have 
made use of the following consideration. In the sum in 
(6), significant contributions come from Sfjiin 
in order to minimize the energy, of which E^^\S) is the 
dominant part. In this range of S, it can be obtained 
that {SaSbSO\S^JSaSbSO) « 6s,s'-iN/4 + 6s,s'+iN/4, 
{SaSbSO\S^\SaSbSO) « -Ss,s'-iN/4 - 5s,s'+iN/4, and 
{SaSbSO\iSf)^\SaSbSO) « Ss', s+2Ny 16 + Ss',sNy 8 + 

In the continuum limit, ipn{S+l)+'4'n{S—'i-) — 2tpn{S) = 
[MS+2) + MS-2)-2i;„{S)]/4: = d''iJn{S)/dS^. There- 
fore one obtains 



(7) 



where u)"^ = 2Kf,/d, f„ = {En + g)/d, with d = {Bb — 
Ba)N/2 - (Ca + Cb + K,- Jz)N^/2, g={Bb- B,)N/2 - 
{Ca + Cb + Ke- Jz)N'^l4.. If > Q, Eq. (7) becomes 
the Schrodinger equation for a simple harmonic oscillator. 



whose n-th eigenvalue is f„ = (n-|- Hence the energy 
level of the original Hamiltonian (4) is given by 



En = {n + l/2)^j2Ked- 



9- 



The interaction gives rise to a term squared in the sum- 
mation variable. This is the reason why the superposition 
coefficients behave like eigenfunctions of harmonic oscilla- 
tors. 

Therefore the energy gap of the Hamiltonian (4) is 



A = ^j2K^d. 



(8) 



Thus 



The ground state of Ti, in 5*^ = sector, is thus 



l*o) = ^V'o(5)|5,0), 



where 



MS)=Ae-^'^^' =Ae 



, 7^'^ 



where the normalization constant A = 
^ s/uj/ erfc{s/ujSmin) , with erfc{x) = dx. 
As ipaiS) is Gaussian, Smin term dominates. The 
smaller the energy gap A, the more dominant Smin term. 
Especially, in case Na = Nb, Smin = 0, thus 



TT 



= (^) 



1/4, 



Hence 



V'o(S') ^ S{S), as A^ 0. 



Therefore, in case Na = Nb, when the energy gap A is 
vanishing while Kg oc l/fi remains finite, the ground state 
of the asymmetric Hamiltonian approaches the symmetric 
state jS* = 0, S'^ = 0). One can say that the symmetry and 
the entangled BEG is protected by the vanishing energy 
gap. 

More generally, no matter whether Na = Nb. as far as 
A <^ Ke, the ground state of the system can be well ap- 
proximated as \Smin,0)- Therefore, the entangled BEG 
persists as the ground state of H, despite that the sym- 
metry is broken in the Hamiltonian. 

A key factor leading to this result is that fl remains fi- 
nite. If, on the contrary, one takes thermodynamic limit 
N oo, oo while N/fl remains constant, then 

Kf. — > 0, consequently tpo{S) becomes independent of S, 
even though one takes A ^ after taking 17 ^ oo. Gon- 
sequently the ground state becomes an equal superposi- 
tion of l^, 0) of all possible values of S, hence breaks the 
symmetry. When the symmetry breaking perturbation is 
infinitesimal, the situation becomes spontaneous symme- 
try breaking (SSB), and is related to Lieb-Mattis infinite- 
range model [6] , with each species of atoms in our model 
corresponding to a sublattice in the latter. 
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(Ba-Bb)-[i^n{S-l) + MS+l)] + iCa + Cb-J, + Ke) — [MS-2) + 2MS)+MS + 2)] = [E^-E^°\S)]MS), 
4 Id 



The present case of persistence of syrametry, and thus 
fragmentation and entanglement, in a finite- volume con- 
densate, could be viewed as a converse case of SSB of 
the ground state of an infinite system. Both arc; due to 
the near degeneracy of the ground state and the low ly- 
ing excited states [7]. For SSB, the limit fl ^ oo should 
be taken before A ^ 0. Likewise, for symmetry persis- 
tence discussed here, should be kept finite before taking 

Consistency of the isotropic parameter point 
with the generalized Gross-Pitaevskii equations. 

- The eff'ective parameters K's are dependent on 

the orbital wave functions pier's, as indicated in 
Eq. (3). These wave functions satisiy the four gen- 
eralized Gross-Pitaevskii equations, which are in turn 
derived from the unique ground state \Go) at the 
isotropic parameter point [1]. The equation for a- 
atoms with spin a is {-^'^'^ + Uaa{r) + [2{N - 
l)/3]5^rVa.(r)P + [{N - l)/3]<7i'rVa*(r)P + [{N - 
l)/3]5^:''Vb.(r)P + m + l)/3]5i"/Vfa*(r)P}'/'a.(r) - 
[{N + 2)/6]ge(l)lif{r)(pba{r)(paa{r) = fJ.aa(f'aa{r), whcrc a ^ 
a, ^aa is the corresponding chemical potential ^, each 
g parameter is the part preceding the integral in the 
definition of the corresponding K parameter in Eq. (3), 
i.e. g^;') = 2TTh^il'^}/t,,„ g, = ^nh^^i'"'^ / ^iab, where 

^e"*"^ = Cjiii = C|ttV equation for (/)6cr(r) is in a 

similar form. 

Therefore, in order to attain the isotropic parameter 
point, both the orbital wave functions and the scattering 
lengths need to be constrained in order that the effective 
parameters satisiy the requirements that Bi = Ci = 
and J2 = Kf,. Because the generalized Gross-Pitaevskii 
equations, which govern the orbital wave functions, are in 
turn derived at the isotropic parameter point, a problem 
arises whether these requirements are consistent with the 
generalized Gross-Pitaevskii equations. 

Here we show that such consistency is indeed guaranteed 
under the following conditions, (i) C/it(r) = Uii{v). (ii) 

The intraspecies scattering lengths ^i"CT2<T3<74's, for phys- 
ically allowed channels, are all equal for each species i, 
denoted as ^j. Note that this means = 2^aa = 
in shorthands, (iii) The interspecies scattering lengths 
satisfy the relations = denoted as ^s"''\ and 

^u''^ = denoted as where the subscripts "s" 

and "d" represent "same" and "different", respectively. 

First, by considering that the nonlinear terms of inter- 
action vanish at the boundary of the condensate, it can be 

^In Ref. [1], as a Lagrange multiplier for the normalization of 
4>a<j, l^acr represents N/2 multiplied by the chemical potential 



seen that C/i|(r) = Uii{v) implies jii^ = ^ii- Then, under 
the above three conditions, in considering that as a wave 
function of bosons, (pi^r is real, one obtains, from the dif- 
ference of the generalized Gross-Pitaevskii equations for 
and (pal, 

N-1 2N+1 (Pb^ (f>bl. N + 2 (Pai-, 

{ 5 9s 5 S'dJlT -7—) = ^ 5e(T -I—)-, 

3 3 (pbl (pbl 6 (pal (pal 

where g^/a = 27r?i^^^°^VMab- Similarly, 
JV-1 2N+1 0„| ,Pai. N + 2 (Pbl (pbi, 

( 9s 5 gd)[- — ) = „ ge(T -I—)- 

3 3 (pal (pal 6 (pbl (pbl 

Consequently, we obtain (pai = (pal and (pbi = (pbi , unless 
{N—l)gs = {2N+l)gd. Subsequently, eii = Cij, according 
to Eq. (2). Consequently, = C, = if the scattering 
lengths satisfy the latter two conditions above. Moreover, 

Transition Temperature. - We briefly discuss the 

issue of transition temperature. A crude estimation can 
be made by following the idea of Ashhab and Leggett on 
SU{2) symmetric model of one species of spin- ^ atoms [8]. 
One may consider four independent interpenetrating gases 
of Nia- atoms of species i with spin a, {i = a,b, a =t 
,1). The transition temperature for each gas is Tjo- « 
S.3lh^m{Nii^/V)'^^^. The error due to ignoring spin ex- 
changes is of the order of VN, and can be neglected in this 
estimation. Hence the transition temperature of the en- 
tangled BEC is roughly the minimum among Tj^'s. Below 
this temperature, the noncondensed atoms of each species 
with each spin is w {mT/3.31h^f/'^ 

Experimental implementations. We now consider 
how to implement this model by using the trapped alkali 
atoms. In order to constrain each atom in the Hilbert 
space of only two spin states, scattering to other spin 
states must be suppressed. 

We propose to represent the two pseudospin states as 
the hyperfine states \F = 2, mp = 2) and \F = 1, mp = 1) 
of an alkali atom with nuclear spin / = 3/2 (^Li, ^^Na, 
^^K, and ^^Rb). Similarly, one can also use \F = 
2,mp = —2) and \F = l^mp = —1) when the magnetic 
field is small so that their energy difference is large. Be- 
cause of the large energy splitting between the two hy- 
perfine states [9], the number of atoms with = 2 is 
forbidden to increase in each scattering. Furthermore, 
in an energetically allowed scattering channel, conserva- 
tions of the total F and mp of the two atoms guaran- 
tee that the hyperfine states of the two scattered atoms 
are either unchanged or exchanged, just as in our model. 
Besides, because of its various advantages, e.g. large in- 
terspecies scattering lengths [10], a mixture of isotopes 



p-4 



Feasibility of Experimental Realization of Entangled Bose-Einstein Condensation iJurop/ii/s. Lett. 86 (2009) 60008 



^^Rb and ^'^Rb can also be used. For ^^Rb, / = 5/2, con- 
sequently |F = 3,niF = 3) and \F = 2, nip — 2) play 
similar roles as \F = 2,mF = 2) and |-F = l,mF = 1) of 
atoms with I = 3/2, respectively. Thus these two sets of 
hyperfine states can represent the two pscudospin states 
of the two species in our model. Similarly, one can use 
\F = 3,mF = -3) and \F = 2,mF = -2) of ^^Rb while 
|F = 2,mF = -2) and \F = l,mF = -1) of *^'^Rb. The 
number of '^^Rb cannot be too large in order to avoid col- 
lapse, as its intraspecies scattering length is negative. 

Magnetic trapping implementation is not favored for the 
following reasons. First of all, only one of the two hyper- 
fine states representing the two pseudospin states is low- 
field seeker and can be trapped. Besides, one of the above 
conditions for the consistency of isotropic parameter point, 
namely Ui-\ = Uii, cannot be satisfied, because the mag- 
netic trapping potential is proportional to gF'nip [9], hence 
has different values for the two hyperfine states. Moreover, 
it is difficult to tune the scattering lengths in a magnetic 
trap. 

In contrast, the implementation can be made in an op- 
tical trap, where the trapping potential is based on the 
single-atom energy shift due to red detuning in coupling 
with the laser, and is independent of atomic spins [9], thus 
the two hyperfine states representing the pseudospin states 
can both be trapped. Moreover, Ui-\ = Uii is satisfied. 
Furthermore, the scattering lengths can be tuned by us- 
ing Feshbach resonances. Therefore, the consistency of the 
isotropic parameter point can be achieved. Even though 
the effective parameters deviate from the exact isotropic 
point, the energy gap can still be tuned to vanish, so that 
the groimd state still approaches the entangled BEC. 

Let us summarize some relevant experiences obtained 
from previous experiments on BEC mixtures. In general, 
for BEC in optical traps, atoms are first prepared in a 
magnetic trap. Optical pumping and adiabatic passage or 
radiofrequency sweep are often used in preparing atoms in 
certain spin states by transferring atoms from other spin 
states. These methods can be adopted in preparing our 
system. In the early experiments in magnetic traps, on 
mixtures of \F = 2, m^? = 2) and \F = l,mF = —1) of 
s^Rb [11] and of \F = 2,mF = 1) and \F = l,mF = -1) 
of ^''Rb [12], spin exchange scattering led to atom loss 
from the magnetic trap. BEC with multiple components of 
different spin states were studied in the optically trapped 
23Na with F = 1 [13] and «^Rb with F = 1 [14], where 
spin exchange scattering indeed did not produce atoms 
with F = 2. Spin exchange scattering were also observed 
in ^''Rb with F = 2 in optical traps, where in some cases, 
e.g. the two initial hyperfine states are \F = 2,mF = 
2) and \F = 2,mF = -2), or both \F = 2,mF = 0), 
scattering between F = 2 atoms lead to F = 1 atoms [15]. 
But such F = 2 to = 1 scattering is prohibited by 
the conservation of total m^^ if the initial hyperfine states 
are \F = 2,mF = 2) and \F = 2,mF = 1), or both 
\F = 2,mF = 2), as in our proposed implementations. 

Two experimental systems of BEC mixtures are very 



close to, and might be extended to, the realization of en- 
tangled BEC. One is a mixture of ^^K and *^^Rb, which 
was regarded as the most favorable candidate for realizing 
BEC mixture [16]. Their BEC mixture has been realized 
in a magnetic trap, with both species in \F = 2,mF = 
2) [17], as well as in an optical trap, with both species 
in \F = l,mF = 1) [18]. In preparing the latter exper- 
iment, atoms are transferred to \F — l,mF = 1) from 
I = 2, mi? = 2) by applying a microwave and a radiofre- 
quency sweep. This procedure could also be used for our 
purpose, with only half of the atoms transferred. Inter- 
species Feshbach resonances, as needed also for our sys- 
tem, were used in this experiment. The other experimental 
system close to ours is that of ^'"'Rb in \F = 2, tuf = —2) 
and ®^Rb in \F = l,mF = —1), in which Feshbach reso- 
nances were used in creating heteronuclear molecules [19] 
and in tuning the two-species BEC [20]. To realize entan- 
gled BEC, we need also ^^Rb in \F = 3, = -3) and 
^^Rb in \F = 2, thf = —2), which can be transferred from 
the other spin state of each isotope by optical pumping or 
radiofrequency sweep. 

As an alternative approach, the two pseudospin states 
could also be implemented by using atoms with negligible 
hyperfine coupling [21]. 

To prepare an entangled BEC, one may first prepare 
four condensates of equal numbers of atoms, of the two 
species with the two spin states, as prescribed above. Then 
conservations of the total spin and its ^-component of any 
two scattered atoms dictate that these atoms remain in 
the manifold of these spin states. The interaction between 
these four condensates builds up coherence among them, 
and realizes entangled BEC of the total system. 

Summary. - To summarize, we have examined the 
feasibility of experimental realization of the entangled 
BEC, i.e. BEC occurring in an interspecies entangled 
two-particle state, as the ground state of a model of two 
species of atoms with spins, proposed in Ref. [1]. The 
entangled BEC is a novel type of fragmented BEC. We 
have analytically shown that in a wide parameter regime, 
the entangled BEC persists as the ground state of the 
concerned model, as far as the energy gap tends to van- 
ish. This makes the entangled BEC more accessible in 
experiments. Subsequently, we established the consis- 
tency between the isotropic point of the effective parame- 
ters and the generalized Gross-Pitaevskii equations, which 
govern the orbital wave functions, on which the effec- 
tive parameters depend. A brief estimation of the tran- 
sition temperature was made. Finally, we discussed how 
to experimentally realize this model with the entangled 
BEC as its ground state. We found that it is suitable 
to use an optical trap. Using two species of atoms with 
/ = 3/2, the two pseudospin states of the two species 
can be realized both in terms of |F = 2, mi? = 2) and 
\F = l,m_F = 1), or both in terms of \F = 2,tof = —2) 
and |F=l,mir = — l).In view of previous experiments, a 
most favorable candidate is the mixture of ^-'^K and *^Rb. 
Another candidate is the mixture of ^^Rb and ^^Rb. For 
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^^Rb, I = 5/2, thus the two pseudospin states are repre- 
sented by |F = 3,mF = 3) and |F = 2,mp = 2) when 
\F = 2,mF = 2) and \F = l,mF = 1) are used for ^^Rb, 
or by \F = 3,mF = —3) and \F = 2,mF = —2) when 
|F = 2,mF = —2) and \F = l,mF = —1) are used for 
87Rb. 

I thank Rukuan Wu for discussions. This work 
is supported by National Science Foundation of China 
(Grant No. 10674030) and Shuguang Project (Grant No. 
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